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ABSTRACT 

Context. While the existence of more than 1800 exoplanets have been confirmed, there is evidence of a wide variety of elemental 
chemical composition, that is to say different metallicities and C/N/O/H ratios. Atmospheres with a high C/0 ratio (above 1) are 
expected to contain a high quantity of hydrocarbons, including heavy molecules (with more than two carbon atoms). To correctly 
study these C-rich atmospheres, a chemical scheme adapted to this composition is necessary. 

Aims. We have implemented a chemical scheme that can describe the kinetics of species with up to six carbon atoms (Co-Ce scheme*). 
This chemical scheme has been developed with combustion specialists and validated by experiments that were conducted on a wide 
range of temperatures (300 - 2500 K) and pressures (0.01 - 100 bar). 

Methods. To determine for which type of studies this enhanced chemical scheme is mandatory, we created a grid of 12 models 
to explore different thermal profiles and C/O ratios. For each of them, we compared the chemical composition determined with a 
C 0 -C 2 chemical scheme (species with up to two carbon atoms) and with the Co-Ce scheme. We also computed synthetic spectra 
corresponding to these 12 models. 

Results. We found no difference in the results obtained with the two schemes when photolyses were excluded from the model, 
regardless of the temperature of the atmosphere. In contrast, differences can appear in the upper atmosphere (P>'''1-10 mbar) when 
there is photochemistry. These differences are found for all the tested pressure-temperature profiles if the C/O ratio is above 1. When 
the C/O ratio of the atmosphere is solar, differences are only found at temperatures lower than 1000 K. The differences linked to the 
use of different chemical schemes have no strong influence on the synthetic spectra. However, with this study, we have confirmed 
C 2 H 2 and HCN as possible tracers of warm C-rich atmospheres. 

Conclusions. The use of this new chemical scheme (instead of the C 0 -C 2 ) is mandatory for modelling atmospheres with a high C/O 
ratio and, in particular, for studying the photochemistry in detail. If the focus is on the synthetic spectra, a smaller scheme may be 
sufficient, because it will be faster in terms of computation time. 

Key words, astrochemistry - planets and satellites: atmospheres - planets and satellites: composition 


1. Introduction 

Observing spectra of transiting exoplanets has proven to be an 
excellent way of obtaining information on the chemical com- 


ysed the Spitzer spectra of WASP-12b. Nevertheless, a C/O 
ratio > 1 in the atmosphere of this planet is quite contro¬ 
versial. Several studies have found that the atmosphere of 


position of their atmospheres. To date, only a few chemical C-rich (Crossfield et alj 20121 

Swain et al. 20131: Mandell et al.l 

species (H, HoU. CO 9 , CH 4 , CO, HCJN, K, JNa) have been in 
ferred thanks to this technique (e.g. ICharbonneau et alJ 2002 

1201 3t IStevensonetal. 12014 

Madhusudhan et alJ l2014ah. To 

’ clearly conclude on the composition of this planet, more ob- 
’ servations with high precision are necessary, using for in- 
“ stance the Wide Field Camera 3 instrument on the Hubble 

Desert et al. 20081 iRedfield et alJ 12008; Swain et al.l l2009alb 

Sing et al.ll201ll: Madhusudhan et al. 2014ah. Accordinglv. om 


is to find the elemental composition and the chemical processes 
that lead to the observed abundances, adjusting different pa¬ 
rameters such as the metallicity, the C/O ele mental abundance 
ratio, or the strength of vertical mixing (i.e. iKopparapu et^ 


___ ng ( __ 

2Qj_2[^ Miller-Ricci Kempton^^l.Jj2()l_d iMoses et al.ll2013alf9 : 

Venot et all2014l : Agundez et al.ll2014bl) . 

The existence of c arbon-rich exoplanet s has bee n suggested 
for the first time by iMadhusudhan et aTI (1201 lal) who anal- 


* This chemical network is available for the community at the 
online database KIDA: Kinetic Database for Astrochemistry at 
http://kida.obs.u-bordeaux 1 .fr/models/ 


iMadhusudhanI (1201 2h and Im OSes et a LdaoBi studied the in¬ 
fluence of the C/O ratio on the chemical composition of hot 
Jupiter atmospheres. They found that models with a C/O ratio 
~ 1 agreed with observational spectra of WASP-12b, XO-lb, 
and CoRoT-2b, which indicates that a wide variety of C/O ra¬ 
tios might be indeed possible in exoplanet atmospheres. 

To study the chemical composition of carbon-rich atmo¬ 
spheres, photochemical models must use chemical schemes 
adapted to this carbon enrichment. When the C/O ratio in¬ 
creases, more complex carbon species are produced. Thus, 
carbon-rich atmospheres are expected to contain heavy hydro- 
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carbons with significant abundances that need to be included in 
the chemical scheme used in the modelling. Most of chemical 
schemes (primarily constructed to study atmospheres with 
solar abundances) do not contain su ch complex spec i es and 
can consequently give biased results. iKopparanu et al.l (1201 2h 
studied the atmosphere of WASP-12b with two different C/O 
ratios (solar and twice solar), but their chemical scheme only 
contained 31 molecules, the heaviest hydrocarbon being C 2 H 2 . 
Thus, the abundance of some hydrocarbons may be over- or 
underestimated because heavier molecu les are absent from 
the chemical scheme. On the other hand, iMoses et al.l (1201 3ah 
used a more complex chemical scheme (90 species), which 
extends up to CeHg, which allows describing the more important 
hydrocarbons. 

Nevertheless, the question of the robustness and the complete¬ 
ness of chemical schemes is not addressed in many publications 
that examine photochemical models. The chemical schemes are 
constructed from re action networks that were originally made 
for J ovian planets (iL iang et all 120031 l2004t iLine et ad l2010l 
et al.ll201 lil2013ai) in which, in the best cases (e.g. 
not for Liang et alT 2003 . l2004l) . high-temperature reactions 
have been added and forward reactions have been reversed to 
reproduce the thermoche mical equilibrium in abs e nce of out- 
of-ea u nibrium proc e sses ( Visscher & Mosesll2011t IVenot et ^ 
l 2012 h . IZahnle et al.l (12009 ) modestly describe the method they 
used to construct their chemical scheme: reaction rates were 
taken from the NIST databas^E and they chose by order of 
decreasing priority between reported reaction rates according 
to relevant temperature range, newe s t revi ew, newest experi¬ 
ment, and newest theory. IVenot et ^ ( 1201 2 h presented the hrst 
chemical scheme that was validated as a whole by experiment. 
Their chemical scheme was constructed in collaboration with 
combustion specialists and was experimentally validated on 
a wide range of temperatures (300 - 2500 K) and pressures 
( 0.001 - 100 bar), making this one of the currently most reliable 
chemical schemes. It is able to describe the kinetics of species 
containing up to two carbon atoms. 


The C 0 -C 2 scheme may be insufhcient, for studying atmo¬ 
spheres enriched in heavier hydrocarbons (with more than two 
carbon atoms). That is why we have developed an extended ver¬ 
sion of the chemical scheme of I Venot et al.l (1201 2 [) , with species 
containing up to six carbon atoms. In Sect.|2]we present this new 
chemical scheme and the grid of models we constructed to test it. 
In Sect.[3]we present the results obtained with the Co-Ce scheme 
and compare them with those obtained with the C 0 -C 2 scheme. 
In Sect. |4] we discuss the implications for interpreting transmis¬ 
sion and emission spectra of exoplanet atmospheres and identify 
species that could play the role of tracers for warm C-rich atmo¬ 
spheres. 


2. Model 

2.1. Chemical scheme 

2.1.1. Reaction network 

The Co-Cg r eaction scheme inc ludes the C 0 -C 2 reaction network 
described in IVenot et al.l (1201 2 b . but also reactions of C 3 -C 6 un¬ 
saturated species, such as C 3 H 6 (propene) or C 4 H 6 (butadiene), 
as well as reactions of small aromatic compounds up to CgHio 
(ethylbenzene). The base of reactions includes the following 
sub-mechanisms: 

' http://kinetics.nist.gov/kinetics 


A primary mechanism including reactions of CvHg 
(toluene) and of C7H7 (benzyl), C6H4CH3 (methylphenyl), 
C 6 H 5 CH 2 OO (peroxybenzyl), C 6 HgCH 20 (alcoxy ben¬ 
zyl), HOC 6 H 4 CH 2 O (hydroxyalcoxybenzyl), OC 6 H 4 CH 3 
(cresoxy), CgHgCHOH, and HOC6H4CH2 (hydroxyben- 
zyls) free radicals. This mechanism is described in 
iBounaceur et al.l (l2005h . 

A secondary mechanism involving the reactions of 
CgHgCHO (benzaldehyde), C6H5CH2OOH (benzyl hy¬ 
droperoxide), HOC 6 H 4 CH 3 (cresol), C 6 H 5 CH 2 OH (ben- 
zylalcohol), CgHm (ethylbenzene), CgHg (styrene), and 
C 14 H 14 (bibenzyp. Th is mechanism is also presented in 
IBounaceur et al.l(l2005h . 

A mechanism for the oxidation of CgHg (benzene) 
(iDa Costa et al.l l2003l) . It includes the reactions of CgHg 
and of cCgH 7 (cyclohexadienyl), cCgHg (phenyl), CgH 502 
(phenylperoxy), cCgHgO (phenoxy), OCgH 40 H (hydrox- 
yphenoxy), cCgHg (cyclopentadienyl), cCgHgO (cyclopen- 
tadienoxy), and CC 5 H 4 OH (hydroxycyclopentadienyl) free 
radicals, as well as the reactions of CgH 402 (ortho- 
benzoquinone), cCgHgOH (phenol), cCgHg (cyclopentadi- 
ene), CC 5 H 4 O (cyclopentadienone), cCgHgOH (cyclopenta- 
dienol), and C 4 H 4 O (vinylketene). 

A mechanism for the oxidation of unsaturated C 3 -C 4 
species. It contains reactions involving C 3 H 2 (propadi- 
enylidene), C 3 H 3 (prop-3-ynyle), the two isomers a-C 3 H 4 
(allene) and P-C 3 H 4 (propyne), the three isomers CgHg 
(allyle), l-CgHg (prop-l-en-2-yle), and 2-C3Hg (prop-1- 
en-l-yle), CgHg (propene), cCgHg (cyclopropene), C 4 H 2 
(diacetylene), the two isomers n-C 4 H 3 (but-l-en-3-ynyle) 
and i-C 4 H 3 (but-l-en-3-yn-2-yle), C 4 H 4 (vinylacetylene), 
the hve isomers n-C 4 Hg (1,3-butadienyle), i-C 4 Hg (1,3- 
butadien-2-yle), 13-C4Hg ( but-l-yn-3-yle), 14-C4Hg (but- 
l-yn-4-yle), and 21 -C 4 Hg (but-2-yn-l-yl), the hve iso¬ 
mers 13-C4Hg (1,3-butadiene), 12 -C 4 Hg (1,2- butadiene), 
cC 4 Hg (methyl-cyclopropene), l-C 4 Hg ( 1 -butyne), and 2 - 
C 4 Hg (2-butyne). This sub-mechanism integrates reactions 
involved in the formation of aromatic compounds and 
has been developed an d validated against exper i mental 
data f rom the literature dWestmoreland et al] 198^ iTsangj 
199 It Miller & MeliusI 1992t iLind stedt & Mauricel 1996t 


Hidaka et al.lll99(C Wang & Frenklaclilll997). 


This database has been used with success in modelling of 
premixed hame of butadiene, propyne, allene, and acetylene 
(iFournet et al.lll99^ and in predicting fo rmation of small aro¬ 
matic compounds for Hame of methane (iGueniche et al.l 1^0091: 
iBelmekki et al.ll^02l) . 

To obtain a good chemical overlap between the reactions bases 
C 0 -C 2 and Cg-Cg, we added a detailed kinetic mechanism for the 
oxidation of a propane/n-butane mixture. This sub-mechanism 
includes a comprehensive primary mechanism, where the only 
molecular reactants considered are the initial organic com¬ 
pounds (here propane and n-butane) and oxygen, and a lumped 
secondary mechanism that contains the reactions consuming the 
molecular products of the primary mechanism that do not react 
in the reaction bases C 0 -C 2 or Cg-Cg. This sub-mechanism has 
been used with success in modelling laminar Hame velocity for 
components of natural gas (iDirrenberger et al]l201 ll) . 

In this reaction base, we included pressure-dependent rate 
constants for unimolecular decomposition, recombination, beta- 
scission, and addition al reactions following the formalism pro¬ 
posed by iTroel (Il974l) . as well as third-body efiiciency coefii- 
cients. As for the C 0 -C 2 chemical scheme, all the reactions of 


2 



















































































Venot et al.: New chemical scheme for studying carbon-rich exoplanet atmospheres 



10-10 10-5 iqO io5 io10 io15 iqO io5 iqIO iqIS io20 


integration time (s) integration time (s) 

Fig. 1: Kinetic evolution of different species with different pressure, temperature, and C/O ratios, computed with the OD model 
using the C 0 -C 2 scheme (full line) or the Cq-Co scheme (dashed line). Kinetic evolutions reach (or evolve toward) thermochemical 
equilibrium (dotted line). 


the Co-Ce scheme are reversed, which allows reproducing ther¬ 
mochemical equilibrium. Most of them (called ’’reversible re¬ 
actions”) are reversed through the principle of microscopic re¬ 
versibility (e.g. Keq = kf/kr, with Keq the equilibrium constant, 
kf the reaction rate of the forwa rd reaction, and k r the reaction 
rate of the reverse reaction, see I Venot et al.ll 20 l^ for more de¬ 
tails), but a few of them (called ’’irreversible reactions”) are re¬ 
versed using experimentally measured reaction rates for both the 
forward and the reverse reaction. This ensures a better reproduc¬ 
tion of the out-of-equilibrium experiments. 

Thermochemical data for molecules or radicals were calcu¬ 
lated and stored as 14 polynomial coeffi cients according to the 
CHEMKIN formalism (iKee et al.lll996h. These data wer e calcu¬ 
lated using the software THERGAS dMulleret al.lll995h . which 
is based on the group and bond additivity methods proposed by 

[BiWi3(IT97^ . 

In summary, the Co-Cg reaction base is composed of 240 
reactants that are involved in 4002 chemical reactions: 1991 
reversible reactions and 20 irreversible reactions. This chemi¬ 
cal scheme is available at th e online database KI DA: Kinetic 
Database for Astrochemistr\Fl llWakelam et alj|20l3) . 

2.1.2. Photolysis 

To model the out-of-equilibrium process that is due to photoly¬ 
sis, we added a set of 113 photodissociations to the Co-Cg chem¬ 
ical scheme. We modelled the stell ar irradiation by using the so¬ 
lar spectrum (iThuillier et alJl2004l) in the range [0-900] nm. 

Eor most molecules, no data at high temperature exist, so 


we used absorption cross-sections at ambient temperature 
that we found in the MPI-Mainz UV/VIS Spectral Atlas 
(iKeller-Rudek et all 1201 3h and in the UV/VisH- Spectra Data 
Basqj. More specific references and explanations of the method¬ 
ology of calc ulating cros s -sectio n s and quantum yield s can 
be found in IVenotetal.l (| 2012 |) . HebrardetalJ (207^ . and 
iDobrijevic et al]~( 2014). Eor CO 2 (IVenot et al. 2013) and NH 3 
dVenot etaL lin prenl) . we used our recent measurements at 
500 K. As we present in Sect. 12.31 we used three different ther¬ 
mal profiles. Eor all of them, we used the same cross-sections 
at 500 K and did not try to adjust the cross-section of CO 2 to 
the exact temperature at which photodissociations occurred be¬ 
cause we are mainly interested in the effect of the new chemical 
scheme at different atmospheric temperatures. However, we ver¬ 
ified with the two different models {T 1000 ^ 0.54 and T 1000 (i.i, see 
Sect. 123] and Table [T] for more explanations and the meaning of 
these symbols) that this approximation has no significant effect 
on the computed chemical abundances. The differences are al¬ 
ways lower than 10%, and for CO 2 , even lower than 0.54% for 
solar C/O. When C/0=1.1, the differences are lower than 30% 
and lower than 6.7% for CO 2 . These differences would not be 
visible on the abundances profiles figures plotted with a loga¬ 
rithmic scale. 


2.2. OD models 

Eirst we compared the kinetic evolution of species calculated 
w ith the two chemica l schemes. We used the OD model presented 
in I Venot et al.l (1201 2 h that allows determining the chemical evo- 


^ http://kida.obs.u-bordeaux 1 .fr/models/ 


^ http://www.uv-spectra.de/ 


3 






























































































Venot et al.: New chemical scheme for studying carbon-rich exoplanet atmospheres 


lution of a mixture at a certain pressure and temperature, with 
no photolysis process. We computed several conditions corre¬ 
sponding to the combination of T = 500, 800, 1000, or 1500 K, 
P - 0.5 or 50 bar, and C/O=0.54 or 1.1, with the C 0 -C 2 and 
the Cfl-Ce scheme. Initial conditions correspond to the thermo¬ 
chemical equilibrium of the bottom level at the profile Tsoo, that 
is, P=1000 bar and 7’=1742 K with C/O=0.54 or 1.1. 

For most species and regardless of the P, T, and C/O conditions, 
the kinetic evolution predicted by the two chemical schemes 
is the same and reaches thermochemical equilibrium. This is 
demonstrated for CO 2 in the first panel of Fig. [T] Nevertheless, 
for some hydrocarbon species, we observe a slight difference in 
the kinetic evolution and sometimes also in the steady state, as 
can be seen in Fig. [1] These differences are found for species 
with more than two carbon atoms, which is normal because the 
C 0 -C 2 scheme is not made to study these heavy species. It con¬ 
tains C„H^ species, with n > 2, only to ensure that species with 
two or fewer carbon atoms will have a correct behaviour. The 
abundances and the kinetic evolutions of these heavy species 
found with the C 0 -C 2 scheme must not be trusted. Thus, it is 
expected that we hnd a difference in evolution for some of these 
hydrocarbons. 


2.3. Grid of 1D models 


T o construct the therma l profil es, we used the analytical model 
of iParmentier & Guillolj ( l2014lPl that was calibrated to match nu- 
me rical P-T prohles o f solar -composition clear-sky atmospheres 


by Parmentier et al. ( 201.'^ . We used the co efficients from 


iParmentier et al.l (1201 5h arid the opacities from [Valencia et alJ 
(l2013h . For simplicity, we did not consider the presence of 
TiO, which may cause a thermal inversion in hot atmospheres. 
Moreover, it is still unclear whether there is thermal inver- 
sion in hot Jupiters (and whether this is caused by TiO) (e.g. 
iMadhusudhan et ah 2014bt Pymentier et al.ll20l4l) . As an illus- 
tration. Ibiamond-Lowe et ^ (l2014t) reanalvsed the Spitzer data 
of HD 209458b and found no evidence for a stratosphere, which 
contradicts previous studies claiming that t he atmosphere of this 
exoplanet presented a thermal inversion (iBurrows et alJ 1200^ 
iKmitson et al.ll200'^ . 

To span a range of conditions relevant to known close-in 
giant exoplanets, we selected as baselines of our study three 
thermal prohles with high-altitude atmospheric temperatures of 
500 K, 1000 K, and 1500 K. These temperature-pressure pro¬ 
hles are shown in Fig. |2 To obtain these prohles, we set the 
irradiation temperature, Tirr, to 784, 1522, and 2303 K. We con¬ 


sidered that jjL - Ij V3, where ^ - cosO and 6 the inclination 
of the stellar irradiation with respect to the local vertical direc¬ 
tion. This choice allowed us to obtain dayside average prohles. 
We considered that the planets have a low internal temperature 
(T/„r = 100 K) and a gravity of 25 m.s“^. At the time of the sub¬ 
mission of our study, a new version of the analytical model was 
released. In this version, the user can no longer directly choose 
the irradiation temperature, but can set the equilibrium temper¬ 
ature, TeqQ, as well as a parameter /, which modulates the hux 
received by the planet. The prohles used in this study correspond 
approximately to dayside average prohles (/ = 0.5) of planets 
with an equilibrium temperature for zero albedo of T^qo - 627, 
1122, and 1718 K for the cool, warm, and hot prohles. Since 
the goal of this study is to compare the C 0 -C 2 and the Co-Ce 
chemical schemes in different conditions, we used the same ther¬ 
mal prohles when considering a solar and not solar C/O ratio. 


^ https://www-n.oca.eu/parmentier/nongrey/nongrey.html 


Table 1: Grid of parameters. 


Parameter 

Range of values 

Symbol 

Temperature 

Cool atmosphere (500 K) 

T 500 


Warm atmosphere (1000 K) 

Tiooo 


Hot atmosphere (1500 K) 

Tisoo 

C/O ratio 

0.54 (solar) 

4o.54 


1.1 

4i.i 

Stellar UV flux 

Irradiation 

pQn 


No irradiation 

Foff 


even if the prohles were calculated assuming a solar composi¬ 
tion. The C/O ratio varies between two values: 0.54 (solar) and 
1.1. Finally, we compared the results including or excluding pho¬ 
tolysis processes. The different parameters used in the models 
are summarised in Table[T]with the corresponding symbols. The 
planet-star distance of each planet has been adjusted to match 
their equilibrium temperature. For the cool, warm, and hot pro¬ 
hles, they have been set to (7=0.2196,0.0549, and 0.0244 AU, re¬ 
spectively. Because many uncertainties exist on the vertical mix¬ 
ing acting in exoplanet atmospheres, we used a constant eddy 
diffusion coefficient, /Qj=10* cm^.s-1, and did not try to adjust 
a complex eddy diffusion prohle. This value is si milar to what 
has b e en used or calculated by previous mo d els (i.e. Lewis et aP 
l20l^ IM OSes et al. I I 20 TTI: iLineetaD I 20 TTI: IVenot et al.1 l2013r 
Howe ver, we recall that this v alue might be too high as pointed 
out by IParmentier et akl (l2013h . We did not explore the space of 
possible values for the eddy diffusion coefficient because it was 
beyond the scope of this paper. Nevertheless, with a stronger ver¬ 
tical mixing, quenching would occur in deeper layers, and for the 
Tiood and Tisoo prohles, the CO/C H4 ratio could beco me lower 
than unity in the entire atmosphere d Venot et al.ll2014ll . 


3. Results 

As can be seen in Figs. [3] and |4] the chemical composition of 
the atmospheres varies depending on the temperature and the 



Fig. 2: Vertical prohles of temperature. The cool atmosphere 
(Tsoo) has an isothermal part from ~20 to 10 “^ mbar at 
500 K(full line). The warm atmosphere (Tiooo) has an isother¬ 
mal part from ~2 to 10“^ mbar at 1000 K (dashed line). Finally, 
the hot atmosphere {T 1500 ) has an isothermal part from ~0.3 to 
10“^ mbar at 1500 K (dotted line). 
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C/O ratio. The photodissociations and the chemical scheme used 
also have an influence. In the following, we comment on these 
differences and explain them for some species. Indeed, by in¬ 
specting the reaction rates, we have identified some formation 
pathways that occur only in the Co-Ce chemical scheme and that 
can explain the departures between the two reactions networks. 
For species for which the two chemical schemes predict different 
mixing ratios, we determined from which reactions these com¬ 
pounds were mainly formed. We compared their reaction rates 
in the two chemical schemes to see which reaction was more ef¬ 
ficient in one of the schemes. Then we examined through which 
reaction the reactants of this latter reaction were formed. We 
proceeded iteratively until we identified a complete formation 
pathway. This identification is a preliminary work. Determining 
the detailed production and destruction pathways of species with 
precision requires the use of algorithms such as th e Pathway 
Analysis Program (lLehmannll 2004 Istock et al.ll 2012 h . These al¬ 
gorithms are very powerful, but need to be adapted to chemical 
schemes as large as our Co-Ce chemical network. For the clar¬ 
ity of Figs. | 3 ] and m we chose to plot only species with signif¬ 
icant abundances (> 10^^), but numerous carbon species have 
an abundance between 10 “'° and 10 “^. We have indicated these 
species in Table | 2 ] 


3 . 1 . Chemical equilibrium 


At chemical equilibrium there are some differences between the 
different thermal profiles and C/O ratios (see the dotted lines 
in Fig. B- Regardless of the C/O ratio, the hotter the atmo¬ 
sphere, the higher the CO/CH4 abundance ratio. CH4 becomes 
more abundant than CO at low temperatures, while CO dom¬ 
inates CH4 at high temperatures. Moreover, the abundance of 
H2O decreases when the temperature increases. We remark that 
when the C/O ratio increases to above unity, then the abun¬ 
dances of carbon-bearing species such as CH4, CH3, and HCN 
increase their abundances significantly. These results agree with 
previous studies on the effect of t he C/O ratio on the chem¬ 
ical atmospheric composition (i.e. iMadhusudhan et^ |201 Ibt 
iMadhusudhanI 1201 ^ iMoses et al.l 1201 3 a ). When H2O is very 
abundant, OH behaves similarly to H, but with a smaller amount. 
However, when the abundance of H2O is quite low (^11 for T 1000 
and T 1500 profiles), the amount of OH is limited by that of water, 
and OH behaves like H2O. Indeed, in each case tested here, the 
main reaction that produces OH is H H- H2O —> OH + H2, which 
means that OH is limited by the reactant of the less abundant 
species. 


3 . 2 . Without photodissociation 

The chemical compositions given by the C0-C2 and the Co- 
Ce chemical schemes are identical (considering species that are 
present in both networks) for all thermal profiles and the two 
C/O ratios when photodissociation is omitted. The only depar¬ 
ture between the two schemes concerns some C„>2Hx species. 
As we explained in Sect. 12.21 the C0-C2 contains species with 
three or four carbon atoms to correctly describe the kinetics of 
the species with two carbon atoms, but the kinetic behaviour 
of these heavier compounds must not be trusted. Thus, it is ex¬ 
pected that the Co-Ce scheme predicts a different behaviour for 
these hydrocarbons. Nevertheless, these small departures do not 
affect the abundances of the main species represented in Fig. [ 3 ] 
The quenching level is different depending on the thermal 
profile. Logically, the cooler the profile, the deeper quenching 


occurs in the atmosphere. For Tsoo, quenching occurs between 
10^ and 10"' mbar, depending on the molecule, whereas for T 1500, 
it occurs much higher in the atmosphere, at about 20 mbar. 


3 . 3 . With photodissociation 
3.3.1. Cool (Tsoo) profile 

For Tsoo, there are differences between the two chemical 
schemes, both for the solar C/O ratio and for the carbon-rich 
case. In both cases, important carbon species such as CO and 
C2H2 are destroyed in the upper atmosphere by photolysis. With 
the Co-Ce scheme, these species are destroyed to a larger extent 
than with the C0-C2 scheme, and we found that the ’’missing” 
carbon is placed in cCeHg through the following chemical path¬ 
way, starting with the photodissociation of CO: 


CO-Hhv 
C-HH 2 
CHH-H2 
3 CH 2 -1- C 
C2H-HH 
2 C 2 H 2 

nC4H3 + C 2 H 2 
ICeHs 
cCeHs -H H 
cCeHe + H 
cCeHv -H H 
CO H- 2 C2H2 -H 2 H2 -I- C 


C -H 0(3p) 
CH-hH 
^CH2 h- H 
C 2 H-HH 
C 2 H 2 

nC4H3 -H H 

ICeHs 

cCeHs 

cCeHe 

cCeHv 

cCeHg 

cCeHg + 0(^P). 


Since it involved heavy hydrocarbons, it can only occur for the 
Co-Ce scheme. 


3.3.2. Warm (Tiogo) and hot (Tisoo) profiles 

As seen in FiglU for any of these two P-T profiles there is no dif¬ 
ference between the results given by the two chemical schemes 
when the C/O ratio is solar, but there are some differences in the 
upper atmosphere when the C/O ratio is > 1 . First, we focus on 
C/ 0 = 1 . 1 . 

In the warm and hot profiles, we found that in contrast to 
the cool profile, cCeHg is unimportant. We discovered that this 
heavy species is not abundant when the temperature is high be¬ 
cause of the reaction nC4H3 -1- C2H2 —> ICeHs, which is the 
limiting reaction in the mechanism identified above for the pro¬ 
duction of cCeHg. The equilibrium constant of this reaction 
(represented in Fig. 13 was calculated as a function of tem¬ 
perature. At low temperature, the equilibrium constant is high 
so kf(T) » kr(T), thus, the addition of nC4H3 with C2H2 is 
favoured to form ICeHs (and then cCeHg). At high temperature, 
the equilibrium constant is very low, so the reverse reaction is 
predominant. The destruction of ICeHs is favoured. As a conse¬ 
quence, the concentration of cCgHg is lower. 

The main species that show differences in their abundances 
when using either the C0-C2 or Co-Ce schemes are CH3, CH4, 
H2O, and CO2, which are more abundant when using the Co-Ce 
scheme than the C0-C2 scheme. C2H2 and H are also slightly 
affected and are less abundant for the Co-Ce scheme. This 
occurs between 10 and 10 “^ mbar for Tisoo, and between 0.3 
and 10 “"' mbar for T looo- We present for each of these species the 
mechanisms we identified to explain these departures between 
the two chemical networks. 
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Fig. 3 : Vertical abundance profiles for the cool (up), warm (middle), and hot (bottom) thermal profiles, with the C/O ratio solar (left) 
and C/ 0 = 1.1 (right). The chemical compositions are calculated with the C0-C2 scheme (dashed line) and the Co-Ce scheme (solid 
line), without photodissociation. Dashed and solid lines overlap, so that the dashed lines are obscured. The chemical equilibrium 
(dotted line) is also represented. 


Methyl radical (CH3): ForCH3, we found that acetylene (C2H2) 
is responsible for the higher mixing ratio obtained with the Co- 
Ce scheme. C2H2 forms CH3 through the following chemical 
pathway, involving heavy hydrocarbons that are only present in 
the Co-Ce scheme; 


C2H2 +H C2H3 

C2H3 + H2 ^ C2H4 + H 
C2H4 + H ^ C2H5 
C2H5 +n 2 CH3 

C2H2 H- 2 H + H2 ^ 2 CH3. 


This explains the lower abundance of C2H2 when the Co-Ce 
scheme is used. 
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Fig. 4; Vertical abundance profiles for the cool (up), warm (middle), and hot (bottom) thermal profiles, with the C/O ratio solar (left) 
and C/0=1.1 (right). The chemical compositions are calculated with the C 0 -C 2 scheme (dashed line) and the Co-Ce scheme (solid 
line), with photodissociation. 


Methane (CH4): Methane is mainly formed (at 99 %) via the 
reaction CH3 -H H2 —» CH4 -H H, so its abundance is linked 
to that of CH3. Because CH3 is more abundant for the Co-Ce 
scheme, there is more methane when using this network. 


Water (H 2 O) and carbon dioxide (CO 2 ): For the Tiooo and 
T 1500 profiles, H 2 O and CO 2 behave in the same way. These two 
species are in eq uilibrium according to the reaction scheme pre¬ 
viously found bv iMoses et al.l (1201 1 ) in hot Jupiters: 


H2O + n <— > OH + H2 
CO -H OH «—> CO2 + H 
H2O H- CO «—> CO2 + H. 

We remark that H2O and CO2 have different abundances depend¬ 
ing on the chemical scheme. These departures are found at ~ 
2xl0~^ mbar for Tiooo and around 10“' mbar for Tisoo. In both 
cases, the two species are more abundant for the Co-Ce than for 
the C0-C2 scheme. Identifying why there is more water when us¬ 
ing the Co-Ce scheme is very difhcult because H2O is involved 
in very many reactions. Nevertheless, we found a possible expla- 
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Table 2 : Highest abundances of some hydrocarbon species not plotted in Fig| 4 ] but with an abundance higher than 10 These 
mixing ratios are obtained using the Co-Ce scheme with the six models including photodissociations. The level pressure (in mbar) 
at which these abundances are reached is noted between parenthesis. Hydrocarbons that are not indicated in this table have an 
abundance always lower than 10“'°. 


Species 

Tisoo ^ 0.54 

Tisoo Ci.i 

T 1000 ( 0.54 

Tiooo Ci.i 

Tsoo ( 0.54 

Tsoo ^ 1.1 

CH 

< 10“‘“ 

4.4x10-* (3x10-") 

< io-‘° 

2.0x10-° (9x10-°) 

< 10-'° 

1.2x10-'° (2x10-*) 

'CHj 

1.7x10-^ (10®) 

7.0x10-° (3x10^) 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 


2.4x10-* (10®) 

1.1x10-° (3x10°) 

< 10-'° 

8.3x10-° (8x10-*) 

2.1x10-° (3x10-*) 

4.6x10-° (3x10-*) 

HCO 

2.5x10-* (10®) 

4.3x10-* (10®) 

1.8x10-'° (10®) 

3.4x10-'° (10®) 

< 10-'° 

< 10-'° 

H2CO 

7.9x10-* (10®) 

1.4x10-° (10®) 

5.7x10-° (10®) 

1.1x10-* (10®) 

1.1x10-° (2x10-*) 

2.8x10-° (2x10-*) 

CH2OH 

1.2x10-“’ (10®) 

2.0x10-'° (10®) 

< 10-'° 

< 10-'° 

< 10-'° 


CH3OH 

LexlO-’ (10®) 

2.7x10-'° (10®) 

4.5x10-'° (10®) 

8.9x10-'° (10®) 

1.8x10-'° (10®) 

3.6x10-'° (10®) 

C2H 

< 10-'° 

1.6x10-° (2x10°) 

< 10-'° 

5.0x10-° (9x10-^) 

< 10-'° 

< 10-'° 

CH2CO 

1.9x10-'° (10®) 

6.8x10-'° (10®) 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

CH2CHO 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

6.5x10-'° (2x10-*) 

1.6x10-° (2x10-*) 

C2H3 

1.1x10-° (10®) 

4.7x10-° (10®) 

< 10-'° 

2.1x10-'° (7x10-*) 

< 10-'“ 

< 10-'° 

C2H4 

8.2x10-* (10®) 

3.6x10-° (10® ) 

8.1x10-° (10®) 

7.1x10-* (9x10-*) 

2.0x10-* (5x10-*) 

7.5x10-* (5x10-*) 

C2H5 

8.1x10-'° (10®) 

3.4x10-° (10®) 

< 10-'° 

1.3x10-'° (10®) 

< 10-'° 

3.6x10-'° (2x10-') 

C2H6 

5.0x10-° (10®) 

2.1x10-* (10®) 

4.1x10-° (10®) 

1.6x10-* (10®) 

6.2x10-* (2x10-') 

2.3x10-® (2x10-') 

aC 3 H 4 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

5.5x10-° (3x10-*) 

1.8x10-* (2x10-*) 

PC3H4 

< 10-'° 

so 

0 

0 

7 

0 

X 

< 10-'° 

< 10-'° 

1.1x10-° (4x10-*) 

3.0x10-° (4x10-*) 

C3H2 

< 10-'° 

6.1x10-° (2x10-°) 

< 10-'° 

1.7x10-° (6x10-*) 

5.3x10-'° (2x10-*) 

1.1x10-° (2x10-*) 

C3H3 

< 10-'° 

2.2x10-° (2x10-°) 

< 10-'° 

7.3x10-° (6x10-*) 

2.5x10-° (2x10-*) 

5.1x10-° (2x10-*) 

C2H3CHO 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

3.4x10-'° (1x10-*) 

8.1x10-'° (1x10-*) 

2C3H5 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

2.6x10-'° (8x10-°) 

3.5x10-'° (8x10-°) 

C3H8 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

4.6x10-° (1x10-*) 

2.9x10-* (2x10-') 

C4H2 

< 10-'° 

1.1x10-'° (10°) 

< 10-'° 

1.5x10-° (5x10-*) 

3.6x10-° (2x10-*) 

6.7x10-° (10-*) 

C4H4O 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

1.4x10-'° (1x10-*) 

IC4H6 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

2.0x10-'° (3x10-*) 

6.9x10-'° (2x10-*) 

iC 4 H 8 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

1.4x10-'° (3x10-*) 

C4H10 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

4.4x10-'° (2x10-') 

cCsHs 

< 10-'° 

< 10-'° 

< 10-'° 

3.1x10-'° (3x10-*) 

1.3x10-'° (2x10-*) 

3.6x10-'° (10-*) 

IC6H4 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

2.6x10-* (3x10-*) 

8.4x10-* (3x10-*) 

cCeHe 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

1.1x10-° (2x10-*) 

3.1x10-° (2x10-*) 

cCfiHsOH 

< 10-'° 

< 10-'° 

< 10-'° 

< 10-'° 

2.2x10-* (3x10-*) 

7.6x10-* (3x10-*) 



Temperature {K) 

Fig. 5 : Equilibrium constant of the reaction 

nC4H3 -I- C2H2 <—» ICeHg as a function of temperature. 
Because the number of reactants and products is different, the 
equilibrium constant is not dimensionless. 


nation: in both chemical schemes, the main reaction (at ~99%) 
producing water is OH -H H2 —> H2O -H H, but the reaction 
OH -H CH4 —> CH3 -H H2O also participates in the formation 
of water. The contribution of this reaction is two orders of mag¬ 
nitude stronger for the Co-Ce scheme than for C0-C2. Moreover, 
methane is ~ 1000 times more abundant at this level for the Co- 


Ce scheme. This can probably explain the higher mixing ratio of 
H2O obtained for the Co-Ce scheme. As we said at the beginning 
of Sect. | 3 ] this hypothesis needs to be confirmed with a specihc 
algorithm, and this will be the subject of a more detailed study. 
The higher amount of CO2 is then a natural consequence of the 
equilibrium between the two species. 


Acetylene (C2H2) We have shown that C2H2 is an important 
species in atmospheres with a high C/O ratio. This species is 
involved in the formation of CH3 and so is indirectly involved 
in the formation of CH4, H2O, and CO2. We can see in Fig. | 4 ] 
that acetylene is a species with a globally high abundance in the 
atmosphere for the ^11 cases. But if we compare this with the 
^0.54 cases, we realise that this is different when the C/O ratio 
is solar. In these cases, C2H2 is not abundant in the atmosphere. 
Furthermore, C2H2 is also abundant in the Tgoo cases in the up¬ 
per atmosphere because of photochemistry. 

As Fig.|6]shows it, for T > 800 K, there is a difference of several 
orders of magnitude between the abundance of C2H2 predicted 
by thermochemical equilibrium when considering a gas mixture 
with a solar C/O ratio or equal to 1 . 1 . The low abundance of 
C2H2 in the /'0.54 cases explains why there is no departure be¬ 
tween the results obtained for the two chemical schemes. If there 
is no C2H2, the chemical pathways leading to the formation of 
CH3 identihed in the paragraph (CH3) cannot occur. CH3, CH4, 
H2O, and CO2 will remain at the same abundances as were found 
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Fig. 6: Abundance of acetylene at the thermochemical equilib¬ 
rium as a function of temperature. Two cases are represented: 
C/O = solar (black) and C/O = 1.1 (red) 


for the C0-C2 scheme. This indicates the important role of acety¬ 
lene in the output of the chemical network calculations. 

4. Spectra 

To determine whether the differences in abundances obtained us¬ 
ing the two chemical schemes have an impact on the spectra, we 
computed the synthetic spectra corresponding t o the 12 atmo¬ 
spheric compositions using the code described in lAgilndez et al.l 
(l 2 Q 14 ah . For this study, the o pacities of C2H 7, C2ii4, and C2H6 
compiled in HITRAN 2012 (iRothman et al.l I^ 13 h were added 
because these C-species can potentially be important in C-rich 
atmospheres under certain conditions. Nevertheless, in the dif¬ 
ferent models we have tested, C2H4 and C2F16 are not sufficiently 
abundant (always lower than 10“^) to affect the synthetic spec¬ 
tra. Note that the opacities of cCeHg and CH3 are not included 
because these data are not known, but they could have an im¬ 
portant influence on the synthetic spectra and should be added 
as soon as they are available. For the same reason, opacities of 
other C„>2H^ species are not included neither. Nevertheless, in 
view of their low abundances (see Table they would proba¬ 
bly not influence the synthetic spectra. We present in this section 
the spectra in transmission that would be obtained during the pri¬ 
mary transit, and then the spectra corresponding to the secondary 
transit, in emission. 

4 . 1. Transmission spectra 

The transmission synthetic spectra computed for the different 
atmospheres are represented in Fig. | 2 l As an example. Fig. [8] 
shows the relative contributions of the main opacity sources for 
the spectra corresponding to the T 1000 thermal profile. 

For the T^qq profile, all the spectra are very similar because 
the atmospheric composition is also quite similar for the two C/O 
ratio cases. The global form of the spectra is due to H2O and 
CH4. Most of the features are caused by the absorption of these 
species. In addition, we note the contribution of NH3 around 1.5 
and 10.5 pm. There are two strong peaks beyond 10 /tm. The first 
one, at 13.6 pm is due to C2H2, and the second one, very close 
at 13.9 pm, is due to HCN. A smaller peak is visible at 14.9 
pm, which is caused by CO2. Despite their similarity, there are 
small departures between the spectra. The differences between 
the C-rich (green, calculated with the Co-Ce scheme) and C/O 





Fig. 7 : Synthetic transmission spectra corresponding to the three 
thermal profiles T500 (up), Tiooo (middle), and Tisoo (bottom), 
represented by the apparent radius. For each thermal profile, 
spectra have been calculated for the atmospheric composi¬ 
tions found using the C0-C2 and Co-Ce schemes for the cases 
C/O = solar and C/O = 1 . 1 , as labelled in each figure. Main 
molecular features are indicated on each spectra. For the T 1000 
and T 1500 cases, features for the C/O solar spectra are plotted in 
orange and features for the C-rich spectra in black. 


solar (red, calculated with the Co-Ce scheme) spectra are due to 
methane, which is more abundant in the ^11 case. But there are 
also differences between spectra corresponding to the two chem¬ 
ical schemes at 3 and 4 pm, between 7 and 8 pm, and between 
13 and 15 pm. The departure at 4 pm is due to the difference of 
the mixing ratio of CO2 obtained when using the two schemes 
(blue vs. red for the ^11 case and green vs. yellow for the ^0.54 
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Fig. 8 ; Contributions of the main opacity sources to the synthetic transmission spectra corresponding to the thermal profile T moo for 
the cases C/O = solar (left) and C/O =1.1 (right). Spectra have been calculated for the atmospheric compositions found using the 
Cq-Co schemes. 


case), and all the others are due to C 2 H 2 , which is more abun¬ 
dant when using the C 0 -C 2 scheme. This means that this species 
contributes more to the transmission spectra. 

For the Tiood profile, we observe strong differences between 
the spectra corresponding to the C-rich and the C/O solar cases. 
First, the global form of the C/O solar spectra is due to H 2 O, 
with contributions of CO 2 around 4.2 and 15 jum, CO around 
2.3 and 4.6 jum, and CH 4 around 2.2 /jm. There is no difference 
between the spectra corresponding to the results obtained with 
the C 0 -C 2 and the Co-Ce schemes. The two C-rich spectra are 
quite different from the C/O solar spectra. Their form is mainly 
due to the H 2 -H 2 collision-induced absorption, with contribu¬ 
tions of several species: CH 4 around 1.16, 1.38, 1.7, 2.3, 3.3 jum 
and between 6 and 9 /tm; HCN around 3 jum and at 13.9 and 
14.2 jum; CO between 4.3 and 5.2 jum; C 2 H 2 around 3 /jm and 
at 7.4, 7.7, and 13.6 jum; NH 3 around 10.5 /tm. HCN is the ab¬ 
sorber that generates the features from 11 to 16 jum, and H 2 O 
is the main absorber from 16 to 30 /tm. Unlike the ^ 0.54 case, 
the two ^1,1 spectra corresponding to the two chemical schemes 
present small departures. We found these deviations between 3 
and 4 //m, between 7 and 8 //m, and between 13 and 15 //m. Like 
for the Tsoo spectra, these departures are caused by the change 
of the abundance of C 2 H 2 between the two chemical schemes, 
but also due to the lower mixing ratio of methane that we obtain 
when using the Co-Ce scheme. 

For the T 1500 profile, we also observe departures between the 
spectra corresponding to the two C/O ratios. For the C/O solar 
spectra, the contributing species are the same as for the Tmoo 
profile: H 2 O everywhere except around 2.3 and 4.6 jim (CO) 
and around 4.2 //m (CO 2 ). The peak of CO 2 at 15 jim is not visi¬ 
ble because of the very much lower abundance of this species in 
the atmosphere. This area of the spectra is entirely dominated by 
water. We note that once again, there is no difference between 
the spectra corresponding to the two chemical schemes, which 
is expected because there is also no difference in the computed 
chemical abundances. The spectra of the C-rich atmosphere are 
different from the C/O solar spectra. Like for the T 1000 profile, 
their form is due to the H 2 -H 2 collision. There are contributions 
of H 2 O around 1.35 jim and from 18 jum until the end of the com¬ 
puted spectra (30 jim). There are many contribution of C 2 H 2 : 
around 1.5, 2.45, 3, 3.7, 7.4, 7.7, and the very high peak at 13.6 
//m. We see CO features at 1. 6 , 2.3, and 4.6 /tm. The contribution 
of CH 4 is very small compared to the two previous cases (Tsoo 


and T looo)- There is a strong peak at 3.32 /rm and a smaller peak 
around 3.49 /tm. Methane is also responsible for the absorption 
peak at 7.39 jim, which is almost mingled with the two strong 
acetylene peaks. Finally, the last contributing molecule is HCN 
with its peaks around 3, 7, 13.9, and 14.2 jum. HCN is also re¬ 
sponsible for the features from 11 to 18 jim. T hese results agree 
with w hat has been found for WASP-12b bv iKoDDaranu et akl 

(I20I2I) . 

4.2. Emission spectra 

Figure|9]presents the synthetic emission spectra obtained for our 
12 atmospheres. As for the transmission spectra, for the cool at¬ 
mosphere, there is almost no difference between the emission 
spectra of the two different C/O ratios, whereas for the two 
warmer profiles the emission spectra are clearly different de¬ 
pending on the C/O ratio of the atmosphere. For the C/O solar 
cases, the spectra are dominated by the absorption of H 2 O for the 
two warmer cases and by water and methane for the Tsoo profile. 
For the C-rich cases, as we explained for transmission spectra in 
Sect. 14.II the spectra are dominated by the H 2 -H 2 collision, with 
contributions of other species, such as C 2 H 2 , HCN, and CH 4 . 

The analysis of these emission and transmission spectra 
shows us that C 2 H 2 and HCN can be considered as good tracers 
of the C/O ratio if the atmosphere presents a sufficiently high 
temperature (above 800 K), thanks to their peaks at 7.4, 7.7, 
13.6, 13.9, and 14.2 pm, respectively, visible when C/O = 1.1. 
For lower temperatures, acetylene and hydrogen cyanide present 
similar mixing ratios in the two C/O cases that we tested and 
thus cannot give in formation on the C/O ratio of the atmosphere. 
As emphasised by iMadhusudhanI (1201 2l) . these two species ab¬ 
solutely need to be considered in spectra models to interpret and 
understand observations of exoplanets atmospheres. 

5. Conclusions 

We have developed a new chemical scheme for studying C-rich 
atmospheres. This Co-Ce scheme describes the kinetics of 240 
species that contain up to six carbon atoms and are involved 
in 4002 reactions. It is available in th e online database KID A: 
Kinetic Database for Astrochemistry (IWakelam et al.ll^l2h at 
http://kida.obs.u-bordeauxl.fr/models/ We used this chemical 
scheme to study atmospheres with different thermal profiles 
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Fig. 9; Synthetic emission spectra corresponding to the three 
thermal prohles T^qq (up), Tioog (middle), and Tisoo (bottom), 
represented by the brightness temperature (K). For each thermal 
prohle, spectra have been calculated with the atmospheric com¬ 
positions found with the C 0 -C 2 and Co-Ce schemes for the cases 
C/O = solar and C/O = 1.1, as labelled in each hgure. 


and C/O ratios and compared the results obtained with those 
obtained with a s maller chemical scheme (C 0 -C 2 , described in 
I Venot et alJ2012h . When photodissociation is neglected, the two 
chemical schemes yield identical results. This strengthens the 
robustness of the small C 0 -C 2 scheme. It perfectly described the 
kinetics of C 2 species, and even a scheme that contains twice 
more reactions gives the same results. 

The introduction of the photodissociations induces, in some 
cases, differences between the two chemical schemes. With 
the warm and hot prohles (~1000 K and ~1500 K), there are 
differences when the C/O ratio of the atmosphere is high (1.1), 


but not when it is solar. This is mainly due to C 2 H 2 , which is 
much more abundant in a C-rich atmosphere and participates 
in the Co-Ce scheme in formation pathways involving heavy 
hydrocarbons. For the cool prohle (~500 K), there are differ¬ 
ences between the two schemes regardless of the C/O ratio 
of the atmosphere. Here again, the high amount of C 2 H 2 in 
the upper atmosphere is responsible for this. Nevertheless, the 
differences obtained using the two chemical schemes do not 
create differences in the synthetic spectra. The analysis of the 
synthetic spectra corresponding to the 12 different cases has 
highlighted the fact the absorption features around 7 and 14 /rm 
due to C 2 H 2 and HCN can be used as tracers of the C/O ratio 
for atmospheres with temperatures higher than 800 K. 

In conclusion, our advice for using of one or the other chem¬ 
ical scheme is as follows: 

- In warm atmospheres (T >1000K) with solar C/O ratios, the 
C 0 -C 2 scheme is sufficient. 

- In cooler atmospheres (regardless of the C/O ratio) or warm 
atmospheres with a C/O ratio higher than 1, the choice of the 
chemical scheme depends on the goal of the study. If the fo¬ 
cus is on computing synthetic spectra, then the use of the C 0 -C 2 
scheme is a reasonable choice because the computation time will 
be shorter. But if the chemical composition is to be studied in de¬ 
tail and the chemical pathways occurring in the atmosphere are 
to be understood, then using the more complete Co-Ce scheme 
is advised. 

We recall that the mixing ratios calculated here could change 
by several orders of magnitude if absorption cross-sections at 
high temperatures are used for all the absorbing species present 
in the chemical network. We only used hot data for CO 2 
(I Venot et al.ll2Q13h and NH 3 (IVenot et al.llin nre^ . but measure¬ 
ments of other species are much needed. In particular, this study 
shows that C 2 H 2 should be the next target of these experiments. 
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